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I.  SUMMARY 


It  has  been  suggested  that  sonic  booms  caused  a  mass  hatching  failure  of  Sooty 
Terns  on  the  Dry  Tortugas,  Florida  [1].  Notwithstanding  that  the  evidence  was 
circumstantial  this  hypothesis  is  often  quoted  as  proof  that  sonic  booms  harm  wildlife  [2], 
Although  subsequent  empirical  studies  with  sonic  boom-like  impulsive  noise  sources 
have  indicated  that  this  thesis  is  unlikely  [2,3]  to  date  no  complementary  analytical  study 
has  been  performed  in  their  support.  This  study  fills  that  gap  by  providing  a 
mathematical  analysis  of  the  problem  of  avian  eggs  subjected  to  sonic  boom 
overpressures  and  interpreting  response  predictions  in  terms  of  egg  damage  potential. 

The  analysis  is  performed  at  two  levels  of  idealization.  The  first  model  is  that  of  a 
spherical  elastic  shell  with  the  embryo  represented  as  an  inertial  concentric  sphere  and  the 
albumen  as  an  acoustic  fluid  that  completely  fills  the  intervening  volume.  The  second, 
higher  fidelity,  model  accounts  for  the  "egg  shape"  and  allows  for  an  air  sac.  In  both 
cases  the  sonic  boom  is  taken  to  be  an  incident  acoustic  wave  with  the  classic  N-wave 
signature.  The  peak  responses  of  the  egg,  specifically  the  peak  shell  stress,  embryo 
acceleration,  and  reactive  force  at  the  substrate  upon  which  the  egg  rests,  are  computed  as 
a  function  of  the  incident  sonic  boom  overpressure  and  compared  with  damage  criteria 
taken  from  the  literature  (Table  5).  It  is  concluded  that  the  overpressures  from  supersonic 
military  operations  with  existing  aircraft,  are  of  insufficient  level  to  cause  avian  egg 
damage  in  general,  and  the  1970  Sooty  Tern  hatching  failure  in  particular. 

II.  ANALYSIS 

The  basic  geometry  for  our  analysis  is  sketched  in  Fig.  1.  A  supersonic  operation 
produces  a  sonic  boom  that  impinges  on  an  egg  resting  on  an  effectively  rigid  substrate, 
or  baffle.  The  overpressure  is  modeled  as  a  plane  acoustic  wave  with  the  classic  N-wave 
signature.  It  is  incident  on  the  egg  and  substrate  with  a  propagation  vector  having  an 
elevation  angle  measured  from  the  vertical  of  a  =  sin'll  / M) .  The  local  Mach  number 
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M=  u/c  with  U the  effective  flight  speed  and  C  the  local  air  sound  speed.  In  equation 
form 


^  '-2-  W  +  1  +  2  ^  • 


t'  -  t-—  [jcsin(a)  -  .ycos(a)] 


with  the  Fourier  spectrum  (transform)  given  by 

jV(co)  =  —  —  [1  +exp(zo)T>]  +  - -  [l-exp(zo)T)] 

2  0)  TO)2 

Pplane(x,y, a, «)  =  exp [  +  ik(x sin(a)  -ycos(a)]  . 

In  the  above,  we  have  assumed  the  harmonic  time  variation  exp(-/a>?),  with  time  t  and 
radial  frequency  CO.  The  variable  k  =  a >/c  is  the  acoustic  wavenumber  and  x  is  the 
duration  of  the  sonic  boom.  N0  represents  the  peak  overpressure  as  typically  measured, 

that  is  at  ground  level  assuming  pressure  doubling. 

Two  levels  of  idealization  are  employed  for  the  egg.  In  both  cases  the  egg  is 
treated  as  a  linear  system  of  continued  and  lumped  elements.  In  one  case  however,  its 
geometry  is  taken  to  be  spherically  symmetric  and  thus  separable,  allowing  for  closed 
form  response  predictions  of  shell  stresses  and  embryo  accelerations.  This  model,  which 
turns  out  to  be  quite  insightful,  is  presented  below  in  Section  A.  A  more  detailed, 
numerical,  model  that  explores  the  asymmetric  influences  of  the  egg  shape  and  air  sac 
characteristic  of  avian  eggs,  is  presented  in  Section  B. 


A.  Analytical  Model  Of  Spherically  Symmetric  Egg 

This  model  is  sketched  in  Fig.  2.  The  shell  is  treated  as  a  thin,  elastic, 

spherical  shell  of  radius  a  and  constant  thickness  hs,  and  the  embryo  as  a  rigid  inertial 
sphere  of  radius  d  and  mass  M  ,  concentric  with  the  shell.  The  albumen  is  represented 

by  an  acoustic  fluid,  of  mass  density  pa  and  sound  speed  ca,  that  completely  fills  the 
intervening  volume.  See  Fig.  2  for  a  sketch  of  the  model  dimensions  and  coordinate 
system.  Here  for  the  sake  of  simplicity,  we  maintain  axisymmetry  by  considering  only 
a  =  0  or  equivalently  M »  1 ,  i.e.,  high  Mach  number  flights. 

The  rigid  substrate  upon  which  the  egg  rests  is  assumed  to  have  two  effects.  It 
doubles  the  incident  pressure  over  the  egg  surface,  and  it  generates  a  reactive  force  that 
constrains  the  radial  motion  of  the  shell  to  be  zero  at  0  —  0 .  This  force,  m  ) ,  is  taken 
to  be  the  resultant  of  a  uniformly  distributed  pressure  over  the  assumed  circular  contact 
area.  A  =  na  2  A2,  where  A  is  one-half  the  subtended  polar  angle.  The  solution  to  this 
problem  is  described  below.  Response  function  expressions,  /(o>),  are  obtained  in  the 
frequency  domain  and,  in  turn,  response  time  histories  are  computed  by  numerically 
taking  their  inverse  Fourier  transforms, 

oo 

f(t)  =  Tf1  Re [J/(co)exp(-  iutydu] 

o 

The  harmonic  equations  of  motion  of  a  thin,  isotropic,  axisymmetric  spherical 
shell  driven  by  the  total  radial  pressure  applied  to  the  shell  p(d, to) ,  may  be  expressed 
as  [4] 

v(e,<o)  =  t  ^4i>»(cos0) 

n-0  wU 

oo 

W(0,CO)  =  £  WnPn( C°S0)  (3 

n-0 
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where  the  modal  amplitudes  F„and  Wn  are  given  by  the  coupled  equations 
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,  A  =  n(n  +  1) 


(4) 


pn  =  [{In  +  l)/2] ■J'  p(d,(D)Pn(cos6)d(cos6 ) 


where  Es  is  Young’s  modulus  and  vs  Poisson’s  ratio  of  the  shell  material,  its  mass 
density,  and  P„( )  is  the  Legendre  polynomial  of  order  n.  Vibration  damping  is 

provided  by  introducing  the  complex  elastic  modulus  Es=>Es(\-ir\s)  with  r^the 
structural  loss  factor.  The  dependent  variables  v  and  w ,  functions  of  the  polar  angle  0, 
are  the  tangential  and  radial  displacement  components  of  the  shell.  These  equations 
include  both  membrane  and  bending  stresses,  the  latter  represented  by  the  p2  terms. 

The  total  radial  pressure  applied  to  the  shell  consists  of  the  external  surface 
pressure  p5(0,co) ,  the  internal  pressure  exerted  by  the  albumen,  pa(0,co) ,  and  the 
reactive  pressure  from  the  substrate  interface,  pR{Q , w) 
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p(0,o>)  =  ^(0,0))  +pa(Q,u)  +Pr(Q,  c*)) 


(5) 


with 

^(0,0)  =  Fr{o)/Ar  0<  A 

=  0  otherwise 

Ignoring  sound  radiation  from  the  surrounding  air,  the  modal  amplitudes  of  the 
external  pressure  exerted  on  the  sphere  which  is  twice  the  incident  pressure  as  given  in 
Eq.  2,  is  [5]. 

Ip  )  =  _?L_  (2n  -  1^-~  N( g>)  ,  (6) 

**H  (ka)2  h'Jka) 

where  h '  (  )  is  the  derivative  with  respect  to  the  argument  of  the  spherical  Hankel 
function  of  the  first  kind  of  order  n. 

The  albumen  being  modeled  as  an  acoustic  fluid,  the  internal  pressure  is  governed 
by  the  acoustic  wave  equation.  The  modal  component  of  the  force  exerted  on  the  shell  by 
the  internal  fluid  may  be  written  as 

(Pa(r))n  =  -  /G)  Z«(r’0))  Wn  ’  a>r>d  »  (7) 

where  the  modal  impedance  Z„  has  a  solution  of  the  form 

Z„(r,o)  =  a  njn  (hr)  +  bnyn{kr)  ,  (8' 

where  jn(  )  and  yn{  )  are  spherical  Bessel  functions  of  order  n,  of  the  first  and  second 
kind,  respectively,  and  k=u>lca.  The  modal  coefficients  an  and  bn  are  defined  by 
enforcing  continuity  with  the  radial  displacement  of  the  shell  at  r=a  and  with  the  rigid 
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spherical  mass  at  r=d.  The  solution  is 


a 

n 

b 

n 

-  «J 

j‘(kp)  y^kt o  >  ‘p  ft 

(/.WWp^W)  [y;W-tpwSk/*>Yn)\  L  0  . 


where  Yn  is  the  translational  admittance  of  the  mass,  given  as 

Y 

n  1  3  <x>My 

Yn*  i-0  • 


’  (9) 


with  My=  py  —  n  d3  where  py  is  the  yolk  mass  density.  To  obtain  the  modal 
coefficients  of  the  internal  surface  pressure  on  the  shell,  Eq.  7  is  evaluated  at  r=a. 

The  modal  coefficients  for  the  interface  pressure  are  obtained  from  the  integral 
expression  in  Eq.  4  as 

=  (2H  +  1)F«(“)^fP  (cos(A))-cos(A)P„(cos(A))]  .  ( 

**•*  2(n+l)  L 

It  remains  to  solve  for  the  value  of  FR(u)IAR  that  enforces  the  boundary  condition 

w(0  =  0,(o)  =  0  i 

Using  equations  (4),  (7),  (1 1)  and  (12)  of  this  reaction  force  becomes 

^  aiz(P«)B 
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(13) 


The  solution  for  the  shell  displacements  may  now  be  obtained  via  substitution  of  Eqs.  (6), 
(7)  and  (13)  into  equations  (4)  and  (5).  The  associated  shell  stresses  are  completely 
defined  by  these  displacements.  Using  conventional  notation,  these  stresses  are  given  by 
[6] 
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w(0 ,  oj) 


•  (16) 


Finally,  the  vertical  embryo  acceleration  becomes 


E  =-  w2  7 


z,-'(r-d,a)Tmmlrr,mJ 


(17) 


B.  Numerical  Model  With  Complex  Egg  Geometry 

With  this  model,  the  geometry  of  the  egg  shell  is  treated  with  greater 
geometric  fidelity  and  the  air  sac  is  represented.  It  is  sketched  in  Fig.  3.  The  shell  is 
again  modeled  as  thin  and  linearly  elastic,  and  capable  of  both  membrane  and  flexural 
stresses.  However  here  the  embryo  and  air  sac,  as  well  as  the  albumen,  are  modeled  as 
acoustic  media.  The  model  was  solved  numerically,  using  finite  elements.  The  egg  shell 
is  developed  from  bending  and  membrane  plate  elements.  Consistent  with  Ref.  7,  the 
shell  geometry  is  characterized  as  a  composite  ellipsoid  described  bj- 
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,  -D<z'<  0 


r'  =  - 
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,  /  — 
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0  <z'<L-D 


0  <  0'<  2ti 


(18) 


where  r’,  z  \  and  0’  are  cylindrical  coordinates,  L  is  the  length  of  the  axis  of  the  egg,  B 
is  the  largest  diameter  perpendicular  to  the  axis  and  D  is  the  distance  from  the  largest 
diameter  to  the  blunt  side  of  the  egg  along  the  axis  (Fig.  3).  Acoustic  fluid  elements 
represent  the  albumen,  egg  yolk  and  air  sac  [8].  It  was  implemented  using  Cosmic 
Nastran  Finite  Element  Analysis  Software.  Our  finite  element  mesh  is  shown  in  Fig.  4. 

The  egg  is  positioned  on  the  substrate  so  that  the  height  of  its  center  of  mass  is 
minimized.  The  substrate  itself  is  again  taken  to  be  rigid  and  as  before  radiation  loading 
is  ignored.  However,  rather  than  simply  doubling  the  incident  pressure  on  the  egg  surface 
as  was  done  for  the  analytical  model,  here  the  surface  pressure  was  actually  computed  for 
a  plane  wave  incident  on  the  egg  in  the  presence  of  the  rigid  planar  boundary,  using  the 
computer  code  Nashua  [9].  The  elastic  egg  is  constrained  such  that  there  are  zero  shell 
displacements  at  the  single  shell-substrate  interface  node.  We  again  interpret  the  reaction 
force  to  be  a  uniformly  distributed  pressure,  now  acting  over  the  average  of  the  areas  of 
the  elements  adjacent  to  the  constrained  node.  The  response  is  obtained  using  Cosmic 
Nastran's  direct  frequency  solution. 

III.  NUMERICAL  RESULTS 

A.  Model  Input  Parameters 

In  this  section  we  describe  the  input  parameters  for  our  models.  For  a  given 
air  sound  speed,  the  sonic  boom,  represented  as  an  N-wave,  is  fully  defined  by  the  peak 
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pressure  and  duration.  Further,  since  our  models  are  linear,  all  response  functions  may  be 
computed  for  a  unit  peak  pressure  and  scaled  accordingly.  For  the  sonic  boom  duration, 
we  take  0.15s,  a  representative  value  for  U.  S.  Air  Force  fighter  aircraft.  Two  incident 
elevation  angles  are  considered,  a  =0°  and  56°,  the  latter  simulating  an  overflight  at 
roushly  Mach  1.2.  These  parameters  are  summarized  in  Table  1. 


Table  1.  N-Wave  Sonic  Boom  Excitation  Parameters 


Description 

Symbol 

Value(s) 

Air  Speed  of  Sound 

c 

330  m/s 

Sonic  Boom  Amplitude 

N0 

1  Pa 

Sonic  Boom  Duration 

T 

0.15  s 

Incident  Angle 
(Mach  Number) 

a 

0°,  56°  (~,  1.2) 

Published  data  [10]  indicate  that  the  basic  structure  and  dimensions  for  eggs  of  a 
variety  of  bird  species,  including  the  sooty  tern,  are  similar  to  those  of  common  chicken 
eggs,  for  which  there  are  a  number  of  published  sources  of  measured  physical  properties. 
Accordingly,  the  material  property  and  geometric  parameters  that  we  have  chosen  for  our 
study,  and  that  are  summarized  in  Tables  2  and  3,  are  based  on  chicken  egg  data. 
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Table  2.  Material  Properties  for  Egg  Models. 


Description 

Symbol 

Value 

Shell  Elastic  Modulus 

Es 

1.1  10"  Pa 

Shell  Poisson’s  Ratio 

Vs 

0.3 

Shell  Mass  Density 

Ps 

2000  kg  /  m3 

Shell  Loss  Factor 

Tls 

0.01 

Albumen  Mass  Density 

Pa 

1000  kg/  m3 

Albumen  Speed  of  Sound 

ca 

1500  m  /  s 

Yolk  Mass  Density 

py 

1100  kg/ m3 

Yolk 

I 

cy 

1500  m  /  s 

The  elastic  modulus  of  the  shell  is  subject  to  considerable  uncertainty,  with  values  in  the 
literature  ranging  from  lxlO10  to  20x10'°  Pa  [1 1,12].  The  value  chosen,  which  falls 
within  this  range,  is  consistent  with  values  computed  from  ultrasonic  measurements  of  the 
compressional  wave  speed  in  a  chicken  egg  shell  that  were  commissioned  for  this  study 
(cs=6,900  m/s),  and  the  listed  shell  density,  Es  =  p SCS2 .  (See  Appendix  A.)  The  mass 

densities  were  estimated  from  measurements  performed  on  purchased  chicken  eggs.  The 
speed  of  sound  of  the  egg  white  and  that  of  the  yolk  were  taken  from  [13].  We  interpret  a 
chick  embryo  to  have  similar  physical  properties  to  a  yolk.  The  shell  loss  factor,  for 
which  no  data  could  be  found,  was  chosen  somewhat  arbitrarily.  Dimensions  were 
estimated  from  measurements  using  chicken  eggs.  The  shell  thickness  listed,  which  was 
measured  at  the  thinnest  part  of  the  shell,  is  consistent  with  other  published  sources 
[7,12].  The  values  of  L,  B,  and  D  agree  with  those  in  Ref.  7. 


Table  3.  Egg  Dimensions 


Description 

Symbol 

Value 

Thickness  of  Shell 

K 

3.6  mm 

Length  of  Shell 

L 

5.7  cm 

Largest  Breadth  of  Shell 

B 

4.2  cm 

Distance  from  Largest  Breadth  to 
Blunt  End  of  Egg 

D 

2.5  cm 

Radius  of  Spherical  Shell 

a 

2.3  cm 

Radius  of  Yolk/Embryo 

d 

1.5  cm 

Area  of  Contact 

Ar 

0.15  cm2 

For  our  spherical  model,  the  radius  a  was  chosen  by  equating  egg  volumes.  For 
the  finite  element  model,  the  shell  interface  area  over  which  the  substrate  reaction  is 
spread  was  conservatively  chosen  to  be  considerably  less  than  that  measured  with  a 
chicken  egg  resting  on  compacted  sand  (7.7  cm2),  while  at  the  same  time  computationally 
convenient. 

B.  Egg  Failure  Criteria 

Peak  egg  shell  stress  and  peak  embryo  acceleration  are  chosen  as  our 
primary  damage  metrics.  Corresponding  failure  criteria,  that  is  maximum  allowable 
levels,  are  shown  in  Table  4.  As  with  its  elastic  modulus,  there  is  a  good  deal  of  spread  in 
the  published  data  for  the  ultimate  egg  shell  stress.  To  err  on  the  conservative  side,  the 
value  of  1 .2xl07  Pa  shown  in  Table  4  is  the  lowest  reported  in  Ref.  1 1,  where  values 
between  1.2xl07  and  5.6xl07  Pa  are  presented.  The  maximum  embryo  acceleration  listed 
is  taken  from  Ref.  14.  In  this  reference,  experiments  are  described  that  indicate  egg 
yolks  subjected  to  the  listed  acceleration  over  a  1.5  ms  period,  yield  no  damage  on  the 
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cellular  level.  The  duration  is  appropriate,  since  as  will  be  seen,  the  computed  embryo 
acceleration  is  a  damped  oscillatory  motion  with  a  period  on  the  order  of  1.  ms. 

Finally,  as  a  result  of  our  study,  a  third  failure  criterion  is  proposed  in  Table  4.  It 
is  the  maximum  allowable  concentrated  radial  force  that  may  be  applied  to  the  egg. 
During  egg  shell  failure  tests  conducted  as  part  of  this  study  by  the  Department  of  Plastics 
Engineering  at  the  University  of  Massachusetts  at  Lowell  (Appendix  B)  it  was  found  that 
egg  failure  correlated  quite  well  with  the  force  applied  by  the  testing  machine.  It  is 
therefore  proposed  as  a  more  robust  failure  metric  than  ultimate  stress.  The  value  shown 
in  Table  4  should  be  conservative  in  that  during  these  tests,  the  eggs  were  placed  between 
rigid  platens,  thus  minimizing  the  area  of  contact  and  in  turn  maximizing  the  stress. 


Table  4.  Egg  Failure  Criteria 


Description 

Value 

Ultimate  Stress  of  Shell 

1.2  107  Pa 

Maximum  Embryo 
Acceleration 

988  g 

(1.5  ms  duration) 

Force  to  Break  Chicken  Egg 

3.2  N 

C.  Response  Predictions 

In  this  section  we  present  computed  response  predictions  for  both  our 
analytical  and  numerical  models.  We  focus  on  the  maximum  egg  shell  stress,  specifically 
the  flexural  stress  at  the  substrate  interface  ( o0  with  z  =  hJ2  and  0  =  0  in  Eq.  14), 
which  we  find  to  predominate.  The  embryo  acceleration,  S ,  and  the  substrate  reactive 
force  on  the  egg,  FR  are  also  presented.  All  results  are  normalized  to  No ,  twice  the 
value  of  the  free  field  peak  overpressure,  to  simulate  the  practice  of  referring  to  the 
strength  of  a  sonic  boom  as  that  measured  at  ground  level,  as  noted  previously. 
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1. 


Time  histories  computed  for  the  shell  stress,  embryo  acceleration  and 
substrate  reactive  force  are  presented  in  Fig.  5.  All  three  signatures  suggest  the  response 
of  a  damped  oscillator  impulsively  executed  at  t=0  and  t=0.15  s,  the  initiation  and 

termination  times  of  the  N-wave.  (Although  not  readily  apparent  in  Fig.  7,  the 
oscillations  correspond  to  a  frequency  of  1,040  Hz.)  To  pursue  this  notion  further, 
consider  the  single  degree  of  freedom,  lumped  parameter  model  illustrated  in  Fig.  6.  For 
the  mass,  M,  we  take  the  total  mass  of  the  egg.  The  spring  stiffress  is  estimated  as  the 
local  stiffness  of  the  spherical  shell  statically  loaded  by  a  uniform  radial  pressure  over  a 
circular  cap  of  area,  AR«  2ahs[l5] 


Eh2  Ar  Ar  2 

K  =  5  —  0.42-0.116  — — +  0.023  — — 

n  ah  ah. 


(19a) 


and  the  dashpot  constant,  B  is  estimated  from  the  structural  loss  factor  in  the  egg  shell 

B  i  TI  JKM  .  (19b) 


Assuming  a  small  loss  factor,  this  results  in  the  natural  frequency 


= 


For  the  parameters  provided  in  Tables  1  and  2  we  predict  a  natural  frequency  of  980  Hz, 
comparing  favorably  with  the  1,040  Hz  oscillations  in  Fig.  5. 

Continuing  with  this  simplified  model,  and  assuming  that  the  egg  is  small  in  terms 

of  acoustic  wavelength  in  air,  that  is  ka  «1,  the  net  translational  force  applied  to  the 
shell  by  twice  the  incident  wave  becomes 
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(21) 


Ftgg(w)  -  i «  — -  2Pbmm (x  =  0,  y  =  0 ,  u)  , 


and  in  turn 


FJ®  “  (2P^=o,y-o,t)) 


(22) 


c  dt 

O 


where  the  location  (x,  y)  is  arbitrarily  chosen  to  be  (0,  0).  The  solution  to  this  forced 
problem  is  readily  obtained.  The  reaction  force,  specifically  the  force  stored  in  the 
spring,  is  given  by 


0  * 

2na3(A 

r 

f 

c 

\ 

exp 

__ 

\ 

2  J 

sin  (mf)  +  exp 


O)  (t  -  x) 


sin(o)  (t-  t)) 


t o  a/c  « 1,  cox  »1 

O  y 


(23) 


The  associated  maximum  stress  in  a  thin  shell,  subject  to  this  reaction  force  is  [15]. 


(  h  ) 

(  a  \ 
ar 

o 

II 

CD 

■i 

ii 

N 

0.42  -  0.282  In 

2’  j 

hs2 

.  ahr 
SJ 

Ar  <  2ahs 


.  (24) 


The  corresponding  acceleration  of  the  mass  (embryo)  is 
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ECO  - 


2tU<23CjO 
^egg  Co 


exp 


2 


sin(o)f)  +  exp 


/  \ 

2 


sin  (go  (t  -  t)) 


) 


to  a/c«  1,  co  x  »1 

0  '  O 


(25) 


Eqs.  23,  24  and  25  are  plotted  in  Fig.  7  for  comparison  with  their  counterparts  in  Fig.  5. 
The  agreement  is  good.  For  example,  the  difference  between  the  computed  peak  stresses 
is  only  17  %.  This  comparison  is  also  reassuringly  insightful  in  that  the  lumped 
parameter  model  shows  explicitly  that  the  peak  amplitudes  are  insensitive  to  the  exact 
choice  of  the  loss  factor,  as  long  as  it  is  small. 

2.  Numerical  Model 

Before  computing  response  time  histories  with  our  numerical  model, 
a  modal  analysis  was  performed.  The  cross-sectional  shape  of  the  fundamental  mode  is 
sketched  in  Fig.  8.  Its  natural  frequency  is  957  Hz.  It  appears  that  most  of  the  egg  mass 
translates  as  a  rigid  body  with  elastic  deformation  limited  to  the  immediate  vicinity  of  the 
substrate  constraint.  The  air  sac  and  yolk  do  not  seem  to  play  a  significant  role.  These 
results  are  consistent  with  the  lumped  parameter  model  of  the  egg  described  earlier.  To 
compare  with  the  findings  from  our  analytical  model,  we  initially  consider  the  case  with 

a  =  0. 

First  for  diagnostic  purposes,  we  follow  the  approach  taken  in  the  previous  section 
to  represent  the  influence  of  the  substrate.  That  is,  we  constrain  the  motion  of  the  egg  at 
the  interface  and  double  the  incident  pressure.  These  results  are  shown  in  Fig.  9.  They 
are  quite  similar  to  our  earlier  predictions.  The  maximum  shell  stress  is  slightly  larger 
than  that  in  Fig.  5.  This  may  be  attributed  to  the  larger  radius  of  curvature  at  the 
interface.  However,  this  picture  changes  with  our  complete  numerical  model.  Here, 
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rather  than  a  priori  doubling  the  incident  pressure  we  account  folly  for  the  diffractive  and 
reflective  effects  of  the  substrate.  These  results  are  shown  in  Fig.  10,  again  with  a  =  0. 
We  note  that  now  the  peak  levels  are  significantly  lower  and  the  time  signatures  more 
closely  resemble  the  incident  N-wave. 

For  insight  into  this  outcome,  we  return  to  the  simple  lumped  parameter  model. 

The  oscillator  was  driven  by  a  net  translational  force  produced  by  the  n=l  modal 
component  of  the  incident  pressure  field  evaluated  on  the  shell  surface.  An  assumed 
doubling  of  the  pressure  doubled  the  force.  However,  with  normal  incidence,  the 
reflected  wave  from  a  rigid  substrate  propagates  in  the  opposite  direction.  Consequently, 
it  tends  to  cancel  any  net  translation  force,  especially  for  small  values  of  ka.  In  turn, 
response  levels  are  reduced.  The  remaining  modal  components  tend  to  produce  local 
reaction  or  spring-like  response,  which  mirrors  the  time  signature  of  the  incident  wave. 
Finally,  in  Fig  1 1  we  present  results  for  the  egg  response  to  an  incident  sonic  boom  with 
a  =5 6°  and  note  that  predicted  levels  are  roughly  50%  lower  than  those  in  Fig.  10. 

D.  Conclusions 

The  ramifications  of  our  predictions  on  failure  potential  are  summarized  in 
Table  5.  Specifically,  we  show  in  Table  5  the  values  of  No  that  would  be  required  for 
the  predicted  failure  metrics  (peak  shell  stress,  egg  acceleration  and  reaction  force)  to 
equal  the  failure  values  in  Table  4.  They  are  unreasonably  high,  especially  in  view  of  our 
consistent  employment  of  conservative  assumptions.  For  perspective,  the  "loudest"  sonic 
boom  recorded  during  a  recent  six  month  military  exercise  involving  385  supersonic 
sorties  was  19.4  psf  [16].  We  therefore  conclude  sonic  boom  overpressures  from 
supersonic  aircraft  operations  are  of  insufficient  levels  to  cause  avian  egg  damage. 


17 


Table  5.  Peak  Sonic  Boom  Overpressure  Required  for  Failure  (a  -  0) 


Failure  Criterion 

N0  Required  For  Failure 

Analytical  Model, 

Pa(psf) 

Numerical  Model 

Pa(psf) 

Maximum  Shell  Stress 

3.75  x  103  (7.8  x  101) 

1.27  x  104  (2.67  x  102) 

Yolk  Acceleration 

4.3  xlO5  (8.98  xlO3) 

3.66  xlO6  (7.65  xlO4) 

Substrate  Force 

2.46  x  103(5.1  x  101) 

1.19  x  104(2.48x  102) 

* 
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Fig.  1  Sonic  boom  model  and  coordinate  system. 
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degrees  of  freedom. 
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Fig-  4  Finite  element  mesh  of  egg. 
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5  Analytical  model  response  to  unit  N-wave  excitation  at 
a=0 .  (a)  reaction  force  (b)  stress  at  position  of  peak 

stress  (c)  embryo  acceleration. 


Fig.  8  Fundamental  mode  shape  of  egg  shell  at  957  Hz. 
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Fig.  9  Finite  element  interim  model  (pressure  doubling  substrate) 
response  to  unit  N-wave  excitation  at  ci=0.  (a)  reaction  „ 

force  (b)  stress  at  position  of  peak  stress  (c)  embryo 
acceleration. 
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Fig.  11  Finite  element  model  response  to  unit  N-wave  excitation 
a=56°.  (a)  reaction  force  (b)  stress  at  position  of  peak 

stress  (c)  embryo  acceleration. 
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APPENDIX  A:  ULTRASOUND  MEASUREMENT  OF 
EGG  SHELL  PROPERTIES 
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INTRASCAN 

Daniel  J.  Cotter 
6  Pleasant  Street 
Easton,  MA  02356 
(508)238-9127 

11/24/96 


Joel  M.  Garreiick 
Kyie  Martini 

Cambridge  Acoustical  Associates,  inc. 
200  Boston  Avenue,  Suite  2500 
Medford,  MA  02155-4243 


Dear  Gentlemen: 

The  acoustic  velocities  of  egg  materials  were  measured  employing  intrascarvs 
Ultrasonic  Materials  Properties  Evaluation  System.  A  Panametrics  M208,  2^,v"_ 
longitudinal  wave  transducer  with  a  delay  line  was  used.  The  velocity  OTSouna  w=^ 
calculated  for  the  egg  white,  yoke,  and  shell  based  on  acquiree  time-oT-nignt  d<*ia. 

Calibration  of  the  Liquid  Path  Distance 

A  simple  measurement  chamber  was  constructed  utilizing  a  common  glass  vial  to 
enableP measurement  of  the  time-of-flight  for  liquids.  The  transducer  was  positioned  to 
a  stop  above  the  liquid  with  the  delay  line  immersed.  Tne  distance  trom  the  end  of  tne 
deiay  line  to  a  reflector  was  measured  employing  the  known  acoustic  ve  ocity  of  ware. 
The  liquid  path  was  determined  to  be  18.9  mm  based  on  the  acoustic  velocity  of  wa.er, 
1.49  mm/microsecond,  and  the  measured  time-oi-flight,  12.7  microsecon  s. 

The  ultrasonic  pulser  receiver  parameters  (gain  time  dependent  gain,  voltage _level, 
etc  )  were  adjusted  to  provide  adequate  signal  to  noise  ratio,  and  allow  ease  o, 
viewing  the  reflected  signals.  The  transducer  was  operated 

Time-of-flight  was  digitally  measured  employing  a  mode  three  technique.  Mode  th.ee 
measures  the  time-of-fiight  between  multiple  back  echoes  to  reduce  oependance  o 
“g and ^mpedanceV.smatch  of  the  transducer.  The  waveforms  were  displayed 
on  a  PC  based  oscilloscope  emulator,  and  time  selective  gates  were  positioned  to 
ensure  discrimination  of  the  reflections  of  interest. 

Determining  the  Acoustic  Velocities  of  the  Egg  White  and  Yoke 

Thp  time-of-fliqht  of  the  egg  white  in  the  liquid  path  of  18.9  mm  was  12.475 
microseconds.  This  corresponds  to  a  longitudinal  wave .velocity _for  the  egg !  «hrt;e  f 
1.520  mm/microsecond.  The  egg  yolk  had  a  time-of-flight  of  12.790  microseconds  o 
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the  same  liquid  path,  or  1.477  mm/microsecond.  The  water,  egg  white,  and  egg 
yolk  were  at  room  temperature  (approximately  25°C)  when  the  measurements  were 
made.  Care  was  taken  to  separate  the  egg  white  and  yolk  before  pouring  the 
materials  into  the  vial.  The  egg  yolk  appeared  to  be  frothy  or  aerated  the 
first  time  it  was  measured.  This  was  not  apparent  on  subsequent  measurements, 
where  the  yolk  did  not  appear  to  "break". 

Measurement  of  the  Acoustic  Velocity  of  the  Egg  Shell 


The  time-of-f light  of  the  longitudinal  wave  though  the  egg  shell  was  measured 
employing  the  same  transducer  previously  described  in  contact  mode  A  liquid 
couplant  was  used  between  the  outside  of  the  egg  shell  and  the  transducer  delay  nne. 
The  time-of-flight  measurements  were  repeatable  and  yielded  59-60  nanoseconds. 
Mechanical  measurements  using  a  micrometer  on  a  flattened  shell  in  proximity  to  t  e 
acoustic  evaluation  site  indicated  406  micrometers  (0.016  inches),  me  calculate 
velocity  was  6.89  mm/microsecond.  This  longitudinal  wave  velocity  is  higher  than 
most  metals  and  more  comparable  to  porcelain  or  ceramic  values  in  the  Lueiature. 


It  was  remarkable  how  repeatable  the  measurements  were  once  the  time  selective 
qates  and  ultrasonic  pulser  parameters  were  established.  The  transducer  could  be 
moved  about  the  egg  shell  while  monitoring  the  digital  signal.  Of  course, jsncettie 
velocity  is  calibrated,  the  thickness  gradient  across  the  egg  could  be  quantnied.  .  he 
Intrascan  system  can  be  operated  as  a  thickness  gage.  The  liquid  and  uquid^leve!  in 
the  egg  could  be  monitored  because  of  the  gross  change  in  impedance  Oi  interlace. 
The  difference  in  reflected  and  transmitted  pressure  could  be  readily  calculated  ior  the 
air  or  liquid  boundaries. 

The  wavelength  of  the  longitudinal  signal  was  on  the  order  of  the  thickness  oMhe 
shell  Ideally  the  specimen  should  be  5  wavelengths  thick.  Higher  irequency 
measurements  could  be  made  for  comparison.  For  example,  Intrascan  owns  and  uses 
a  Panametrics  5218  thickness  gage  pulser  receiver  with  a  100  MHz  bandwidth  for 
evaluation  of  thin  materials.  The  shear  wave  velocity  could  not  be  readily  measured. 

The  return  signal  of  a  normal  incidence  shear  wave  transducer  was  not  detectable,  it 
is  probable  that  the  signal  was  lost  in  the  “main  bang”  (excitation  pulse  or  front  surface 
interface).  The  velocity  of  the  shear  wave  is  typically  50%  of  the  longitudinal  wave; 
however  a  5-MHz  normal  incidence  shear  wave  transducer  was  employed 
Consequently,  the  wavelength  was  approximately  doubled.  The  appropriate  first 
critical  angle  for  the  transducer  can  be  calculated  based  on  the  measured  velocities 
for  the  shell.  Additionally,  higher  frequency  probes  are  commercially  available. 
Measurement  of  the  shear  wave  velocity  and  the  calculation  of  the  elastic  properties 
for  the  egg  shell  were  beyond  the  scope  of  the  initial  effort. 
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A  three  part  proposal  wiii  be  prepared  under  a  separate  cover. 

1)  Assessment  of  the  feasibility  of  measuring  the  shear  wave  velocity, 
thereby  enabling  calculation  of  the  elastic  properties  (Young  s  Modulus, 
Shear  Modulus,  and  Poisson’s  Ratio)  of  the  shell. 

2)  Assessment  of  the  feasibility  of  monitoring  the  impact  of  an  equivalent 
sonic  pressure  pulse  (sonic  boom)  on  the  egg  liquids  in  real-time.  The 
compression  and  rarefication  of  a  liquid  has  been  measured  in  real-iime  in^ 
Intrascan’s  laboratory.  This  is  typically  done  to  characterize  the  wavefront  of 
puises  from  transducers  in  water. 

3)  Quotation  and  description  of  Intrascan’s  system. 


If  you  have  any  questions,  please  call. 


Sincerely, 

/  -  ' 

V  : 

Daniel  J.  Cotter 
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APPENDIX  B:  MECHANICAL  STRENGTH  TEST  RESULTS  FOR 

EGGS  AND  EGGSHELL 


Prepared  by: 

Plastics  Engineering  Department  of  the 
University  of  Massachusetts  at  Lowell 

(Fig.B.l) 


37 


DEC  Its  '  tib  k flfii*.  j>njr-n  r-rcii'-J  ^ 


Cambridge  Acoustical 
Summary  of  Results 

Dec*  14,  1996 
Page  1  of  2 


Rate  of  extension  for  tests  unless  otherwise  noted: 
rate  =  0.2  inches  per  minute 

Eaq  puncture  -  Pressure  foot  area  =  0.312  cm  The  load  xs 
expressed  as  the  maximum  load  to  first  failure  point.  This 
represented  the  point  where  the  shell  started  to  crack. 


Whole  egg  (full)  -  load  to  initial  failure  -  pounds 


7.58 

(3-44) 

7.46 

(3.39). 

7.18 

(3.26) 

Ave.  7.41 

(3.36) 

Blown  (empty)  egg  shell 

9.02 

(4.10) 

7.90 

(3.59) 

3.88 

(1.76) 

8.82 

(4.00) 

6.81 

(3.09) 

8.85 

(4.02). 

5,30 

(2.41) 

Ave.  7.23 

(3.28) 

Slope  of  load/ compression  curve: 
"elastic  region"  =  1,158  lbs/inch 
at  br?ak  point  =■  1,371  lbs/  inch 
load  rata  -  0.02  IPM 


-  load  to  failure  -  pounds  (kg) 


Slope  of  load/compression  curve: 
"elastic  region"  =  769  lbs/ inch 
this  value  near  constant  over  full 
test  range.  load  rate  -  0.02  IPM 


Eaa  crush  -  The  load  is  expressed  as  the  maximum  load  exerted 

between  two  platens,  platans  larger  in  /Shere  thS9lhell 

first  failure  point.  This  represented  tne  point  where  the  sneii 

started  to  crack. 

Whole  egg  (full)  -  load  to  initial  failure  -  pounds  (kg) 


6.26 

(2.84) 

7.34 

(3.33) 

6.67 

(3.03) 

Ave.  6.76 

(3.07) 

B.l 
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Cambridge  Acoustical 
Summary  of  Test  Results 
Dec.  14,  1996 
Page  2  of  2 

Hoop  stress  -  The  egg  shell  was  sliced  to  form  a  ring  aroiind  the 
largest  diameter  of  the  egg.  The  ring  was  tested  on  a  fixture 
under  tensile  load.  The  fixture  was  "fitted"  into  the  shell  ring 
to  restrict  deformation  of  the  shell  during  extension. 


Hoop  width  ( in) : 


0.477 

0.473 

0.472 


0.472 

0.473 

0.474 


0.476 

0.476 


0.472 

0.476 


Ave.  =  0.475  inches  (1.207  cm) 

Hoop  Stress  -  pounds  (kg)  Actual  load  is  one  half  the 
reported  values  listed  below. 

1.24  (0.565) 

1.41  (0.640) 

2.23  (1.035) 

1.93  (0.900) 

1.62  (0.734) 

Ave. 1.71  (0.775) 

Hoop  compression  -  The  egg  shell  was  sliced  to  form  a  ring  around 
»h/i ^raest-diameter  g£_the  egg.  The  ring  was  tested  between  two 
flat  platens  under  a  compression  load.  The  egg  shell  was  placed  cn 
a  small  dab  of  plaster  of  paris  to  prevent  the  ring  from  slipping 
during  test.  The  shell  ring  to  free  to  deform  during  compression. 

Eoop  width  (in):  0.471  0.476 

0.472  0.477  0.472  0.473 

0.476  0.474  0.476  0.472 

n  a?*  0.475  0.474  0.475 


0.471 

0.472 

0.476 

0.476 

0.474 


0.476 

0.477 

0.474 

0.475 

0.477 


0.476 

0.472 

0.476 

0.474 


Ave.  =  0.474  inches  (1.204  cm) 


Hocp  compression  -  pounds  (kg) 


-  "Elastic"  slope  lb/ inch 


#T 

0.339 

(0.154) 

0.754 

#1 

0.293 

(0.133) 

0.692 

#2 

0.229 

(0.104) 

0.416 

#3 

0.326 

(0.148) 

0.753 

#4 

0.282 

(0.128) 

0.661 

#5 

0.313 

(0.142) 

0.759 

Ave. 

0.297 

(0-135) 

0.673 
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Fig.  B.l  Egg  strength  testing 
apparatus . 


(a)  Sample  Configuration  For 
Measuring  Egg  Failure  Load 


(b)  Sample 

Configuration 
For  Measuring 
Egg  Shell 
Failure  Stress 


